High-fidelity ion-trap quantum computing with hyperfine clock states 
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We propose the implementation of a geometric-phase gate on magnetic-field-insensitive qubits with a z - 
dependent forces for trapped ion quantum computing. The force is exerted by two laser beams in a Raman 
configuration. Qubit- state dependency is achieved by a small frequency detuning from the virtually-excited 
state. Ion species with excited states of long radiative lifetimes are used to reduce the chance of a spontaneous 
photon emission to less than 10 -8 per gate-run. This eliminates the main source of gate infidelity of previous 
implementations. With this scheme it seems possible to reach the fault tolerant threshold. 



PACS numbers: 03.67.Lx, 03.67.Pp, 32.80.Qk 

Introduction. One of the big challenges in quantum com- 
puting today is to perform nearly perfect gate operations. With 
trapped ions, initialization, read-out and single qubit opera- 
tions have been already demonstrated by various groups with 
very high fidelity. The implementation of a two-qubit gate 
operation is more involved as it demands coupling the inter- 
nal states of two ions separated by a distance many orders 
of magnitude larger than the Bohr radius. Ideally, the gate 
should operate at high speed and with a fidelity allowing for 
fault- tolerant quantum computation. To comply with the latter 
requirement, the gate operation should not couple the qubits 
dissipatively to the environment (for example, by spontaneous 
emission), nor should the qubits be affected by dephasing 
mechanisms like fluctuating magnetic fields or path lengths 
in the optical setup. With trapped ions impressive progress 
towards high-fidelity two-qubit gates has been made over the 
last couple of years, yet all current experimental realizations 
IlllfiSH] achieve fidelities not exceeding 97% and thus still 
fall short of coming close to the desired level of precision 0]. 

Optical Raman fields are a convenient tool to create strong 
spin-dependent forces on hyperfine ion-qubits OL S & 0], 
where "spin" refers to the effective Pauli spin associated with 
the qubit 's two level system. Acting on a single ion, these 
forces can entangle the internal- spin and motional degrees 
of freedom of the ion |6l [ZD ; following disentanglement, the 
ion acquires a net geometric phase that is spin-dependent (H]. 
When the force is exerted simultaneously on two ions, the 
geometric phase depends nonlinearly on their spins and can 
be used to entangle the ions' internal degrees of freedom 
OL HI IE HI OH- Quantum gates based on these collective 
laser-ion interactions come in two flavours. In the first ap- 
proach, experimentally pursued by Leibfried et. ah flJJ], a 
moving standing wave pattern of two laser beams induces, for 
properly chosen laser beam polarizations, a spatially-varying 
state-dependent AC-Stark shift. This generates in turn a (in- 
dependent force, i. e. , a force on the ions whose amplitude de- 
pends on the internal energy states. The resulting controlled- 
phase gate, called a 2 -gate, is remarkably fast and robust. 

As the Raman coupling in [1] is mediated by a short-lived 



Pi/2 level (see Fig. Q]a)), spontaneous emission dominates 
the error budget even for a large detuning A from the virtual 
level full . Unfortunately, the use of a large Raman detun- 
ning makes the gate operation very inefficient when acting on 
"clock" states (states whose energy splitting is first-order in- 
dependent of changes of the magnetic field) [Q, E3, EH], for 
no differential Stark shift can be induced on them in the limit 
A ^> ujq 1 13] , where c^o is the hyperfine structure splitting 
of the qubit states. Therefore, the remarkably long coherence 
times 1 14] available for clock states cannot directly be com- 
bined with the intrinsic robustness of the a 2 -gate when based 
on Raman transitions with electric dipole coupling. 

The second type of collective quantum gates, pioneered 
by S0rensen and M0lmer @], entangles the qubits by induc- 
ing collective spin flips of both ions. It is based on a (in- 
dependent force, with being a linear combination of the 
Pauli spin operators a x and a y , capable of operating on clock 
states [7]. Even though this gate is formally equivalent to 
the a 2 -gate in a rotated basis though, it has the disadvantage 
that it involves spin-flips, which make the gate more sensi- 
tive to magnetic field fluctuations than the a 2 -gate Hi. Spin 
flips also constitutes a drawback for quantum computation in 
decoherence-free subspaces protected against collective de- 
phasin g, fo r which controlled-phase gates are naturally better 
suited 11511 . Furthermore, the required laser beam configura- 
tion gives rise to strong sensitivity of the qubits' coherences 
on the optical phases of the driving fields, constituting a seri- 
ous limiting factor to the fidelity of the gate fl6ll unless spe- 
cial spin echo techniques are applied within the gate operation 
10, 0. In addition, for this type of gates spontaneous emis- 
sion also sets a fundamental major limit to the fidelity. 

In this paper, we investigate how to overcome the a z - 
gate's spontaneous emission problem and its inefficiency with 
clock states by using electric quadrupole transitions for me- 
diating the Raman coupling. The main idea is to induce 
a state-dependent force on the ions via Raman lasers tuned 
close to a narrow transition. We illustrate the idea with the 
Si/ 2 <-> D5/2 transition in 43 Ca + , but it can also be ap- 
plied to other ion species with a similar level structure fl7ll . 
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Two non-copropagating lasers are tuned such that the two 
qubit levels, | |) = \S 1/2 (F = 3,ra F = 0)) and | |) = 
I S1/2 (F = 4, rnp = 0)), experience opposite detunings from 
the metastable D 5 / 2 l eve l ( see Fig- HI &))• Two ions placed 
inside the standing wave pattern created by such laser beams 
experience a force. The direction and amplitude of the force 
on each ion depends jointly on the ions' internal state con- 
figuration, the distance between them, and the collective mo- 
tional mode chosen. As in Ref. yj], a conditional dynamics is 
achieved by introducing a frequency difference of both beams 
close to the frequency a vibrational mode of the ion string. As 
a result, the state-dependent force becomes periodic and ex- 
cites, for instance, the center-of-mass (CM) mode if and only 
if the two ions are in opposite internal states and separated 
by a semi-integer multiple of wavelengths. The periodicity 
induced by the detuning from the vibrational frequency en- 
sures that the ion string eventually returns to its initial mo- 
tional state. During this process, however, the ion string picks 
up a phase proportional to its excitiation, which allows to re- 
alize a maximally entangling controlled-phase gate. 

Distinguishing between the internal ionic levels only via 
their respective detunings works for any pair of qubit levels 
with a sufficiently large energy difference regardless of their 
magnetic properties, in particular for clock states. However, it 
requires that the mediator level (the D 5 / 2 l eve l m our exam- 
ple) is long-lived, to achieve a reasonably small spontaneous 
emission rate. Concentrating again on 43 Ca + , the D$/ 2 level 
has a lifetime of about 1 s, such that the probability of sponta- 
neous emission of a photon from both ions, and for a Raman 
detuning of half the energy splitting of the groundstate hy- 
perfine structure, is well below the asymptotic threshold for 
fault- tolerant quantum computation found in . 

Spectral discrimination. We now describe the emergence 
of the state dependent force via spectral discrimination and 
the implementation of the a 2 -gate. For a more detailed expla- 
nation on general spin-dependent forces we refer the reader 
to Ref. II13I1 and references therein. We are interested in the 
implementation of the interaction Hamiltonian 

+ h.c.)\rrii)(mi\, (1) 
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with mi =| or I being the pseudo-spin state of ion z, for 
i = 1 or 2. Hamiltonian (Q]) describes a particle with internal 
states I mi) trapped in a harmonic potencial of frequency v cm 
and driven by a classical periodic force of amplitude Fo m . 
and frequency 6 <C v C m- We consider the collective center- 
of-mass (CM) motional mode as the harmonic oscillator be- 
cause it exhibits lower levels of dephasing [18]. Nevertheless, 
the treatment can be straightforwardly applied to the stretch 
mode, too. In Eq. (Q]), z Ucm = ^/h/2Mv cm , where M is the 
single-ion mass, is the root mean square spread of the mo- 
tional ground state wave packet and aj m is the creation oper- 
ator of one phonon of the CM mode. Upon evolution with 
Hamiltonian (QJ the oscillator is periodically driven through 
a circular trajectory in phase space. For an evolution time 
T = 2nirS~ 1 , with n G N, it ends up back in its initial posi- 
tion, but with each state |ra$) having acquired a geometrical 

phase $(T = 2n7rS- 1 ) = 2nix\ Z -^^ 
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Figure 1: (Not to scale, color online). The usual implementation a) 
uses a P1/2 level as the excited state and works with A ^> cjo; the 
qubit states are magnetic-field- sensitive and the differential coupling 
is achieved mainly through polarization discrimination. In the new 
implementation b) the excited state is a D 5 / 2 metastable level and 
A is small enough for the qubit states to be spectrally discriminated, 
which allows for the use of clock states. Maximal discrimination 
occurs at A = where the effective coupling of | f) to its motional 
sideband is exactly opposite to that of | j). 



is exerted simultaneously on both ions, the phase eventually 
acquired by each composite state | mi , 777-2) can be set to yield 
the maximally-entangling a 2 -gate by properly controlling the 
values of Fq 19 i*b m2)2 and S, as is described below. 

Fig. [T] b) sketches the relevant frequencies, 
00 a = ^t,e + A — (z/ cm — S) and 00 b = ^t, e + A, 
uo b — = v cm — S, with oo^ e being the frequency difference 
between the excited state \e) and | |), and c^o between | j) 
and ||). We have A ^> is cm ^> S, while A and c^o are of 
the same order of magnitude. Since the laser detuning A is 
still much larger than the D 5 / 2 state's hyperfine structure we 
can treat the latter as an effective single-level without internal 
structure. The individual interaction Hamiltonian in the 
Schrodinger picture and optical rotating-wave approximation 
(RWA) are then given by H[ = ^Ez,m, 9i,i,rnM){mi\ 
gilfci.^CaL+^-^t-w+fci,^] + Hc j]^ xhe Rabi 

frequency gi y i ymi , of laser I = A or B, couples level 
1 777,$) with |e) of ion i, k\ z is the component of the /-th 
laser's wave vector along the trap axis z, ipi its optical 
phase, and zo i the equilibrium position of ion i. If (I) 
\gi,i^\ < A and \gi^i\ <C |A-cj |, and since A w \ A-oo \ 
is much larger than the excited state's linewidth, level 
|e) can be adiabatically eliminated. The latter, together 
with the RWA eliminating terms rotating at frequencies 
A and A — ooq yields the following effective Hamilto- 
nians in the interaction picture: H[ = h^ m . [Xrrn,i + 
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time-dependent components. 

The last Hamiltonian still contains fast oscillations at fre- 
quency v cm . In (II) the resolved side-band limit \0 mi ^\ <C v cm 
only the stationary (non-oscillating) terms give a significant 
contribution to the effective Hamiltonian, whereas the other 
can be neglected under the RWA. Furthermore, (III) in the 
Lamb-Dicke limit (LDL) ?/ c 2 m (n cm + 1/2) < 1, where n cm is 
the mean population of the CM mode, only the terms linear 
in r] cm contribute. Then the effective Hamiltonians are finally 
given by: 

HinU = [Xmi,i 
mi 

+ {iO m ^Vcrne-^at m e* 6t + h.c)] |m;)(m;|. (2) 

Let us momentarily disregard the time-averaged Stark shift 
from Hamiltonian (|2]) and consider only its time-dependent 
component that contains the vibrational operators generat- 
ing the interaction between both qubits. Under this assump- 
tion it suffices to define F 0mi % z Vcm = ihO mii ir] cm e~ lcf)i to 
identify Hamiltonian © with (Q]), realizing thus the desired 
a z -dependent force on both ions. For simplicity, we take 

9l,l,m = 9l,2,m = gi,m => @rn,l = #m,2 = #m |^0 m ,i | = 

l^b m , 2 1 — \Fo m | , as no individual laser addressing has been 
assumed. The extension to non-equal coupling though is 
straightforward. The phase difference between the forces ap- 
plied to both ions is then determined by the ion spacing at 
equilibrium: <pi — <p2 = Ak z Azo, with Azo = zq 1 — zo 2 . 
Here it is convenient to set the spacing equal to a semi-integer 
multiple of the optical wavelength: Ak z Azq = (2n + l)7r, so 
that both ions experience opposite forces. In this way the total 
force on the CM mode F 0mi ^ 2 z Ucm (F 0mi + F 0m2 )z Ucm = 
)Vcme 1(1)1 is non- vanishing only for anti-aligned 

spins. 

After one phase-space round trip of duration T = 2ttS~ 1 
(the shortest gate-time possible) the composite states evolve 
as: |rai,77i2) — > |mi, 7712), for mi = rri2', and |mi, UI2) — > 
e^|mi, m 2 ), for mi 7^ m 2 , with $ = 2ir\(6^ - ^m^ -1 1 2 . 
Maximal entanglement is achieved for = | ; and the bigger 
the difference between 0^ and #| , the faster the gate. Imposing 
thus <1> = 77 and using the definition of m , the general dis- 
crimination condition for maximal entanglement is obtained: 



(3) 



Since we are most interested in working with qubit states with 
the same magnetic properties no discrimination via polariza- 
tion can take place, i. e., g\^ = gi^ = gu for I = A ox 
B. Therefore, the only way to fulfill condition ©, which 
then reads \gBg\ [A -1 — (A — ^o)" 1 ] | = | 2^ I' is mrou gh 
spectral discrimination. Maximal discrimination occurs at 
A = when 0^ = —6^. As shown in figure [T] b), for this 
detuning, | |) couples to the excited state \e) through blue- 
detuned light and | j) through red-detuned light of the same 
intensity. So both states' effective couplings with the motional 
sideband (and therefore also the forces each spin state experi- 
ences) are exactly opposite. Thus the optimal form of condi- 
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For 43 Ca + , the ground state hyperfine splitting is ujq = 
2tt x 3.226 GHz. The resulting Raman detuning is 20 times 
larger than the hyperfine splitting of the D 5 / 2 manifold, so 
that the single-level approximation is comfortably satisfied. It 
is also much larger than the linewidth 7^ = 2ir x 0. 1 8 Hz. For 
the typical experimental parameters: v cm = 2tt x 1.2 MHz, 
Vcm = 0.1 and 5 = 2tt x 1 kHz, and taking \g^\ = = \g\, 
the value g = 2tt x 2 MHz is obtained from Eq. © for the opti- 
cal quadrupole couplings. All the other approximations made 
in the derivation of Hamiltonian © are in turn well satisfied. 
It is: (I) \g\ = 27rx2MHz<27rxl.613GHz= |A-cj | = 
A; (II) = 0.002 < 1; and (III) r/ c 2 m (n cm + 1/2) = 
0.01 (n cm + 1/2), which is much smaller than 1 for the mode 
populations of interest. On the other hand, the total duration 
of the gate (limited by the available amount of laser power) is 
T =1 ms, which is short as compared to the long coherence 
times expected for clock states 111411 . To speed up the gate time 
to 100 [is for example, a total laser power of about 1 W, when 
focussing the laser beams to waist sizes of 6 fim, would be 
required. 

Not all types of spontaneous emission affect the gate's per- 
formance equally [20]. The dominant source of gate error 
comes from inelastic Raman scattering that affects the ion- 
qubit coherences directly, whereas elastic Rayleigh scattering 
can only affect the gate's fidelity indirectly through the ions' 
recoils. In what follows, we calculate the total probability of 
spontaneous photon scattering and simply take it as an overes- 
timated figure of merit for the gate error. The mean probability 
of off-resonant direct excitation of the excited state \D 5 / 2 ) un- 
der the action of the two ^--detuned lasers is approximately 

Poff = = 2 x 10 -6 . The total probability of spontaneous 

emission from both ions during the gate-time T = 2n5~ 1 is 

then: p T = 2 x PoS x lD x 2^ = 2 x 8^ x lD x = 

izLDLo = q 3 x 10 -9 , where condition © has been used. 

?7cm ^0 

This value is more than four orders of magnitude below the 
asymptotic threshold value set in [5] for fault-tolerant quan- 
tum computation. When it comes to choosing a specific ion 
as a quantum information carrier in the case of gate opera- 
tions mediated by a short-lived P level, it seems favourable 
to avoid ion species with low-lying D levels [2(|. The gate 
scheme proposed here not only works for any ion species with 
a low-lying narrow transition, but also it is exactly the use of 
the low-lying narrow transition as the mediator what enables 
an error threshold four orders of magnitude smaller than for 
gates mediated by short-lived P levels. 

Single -ion phase compensation. Let us now go back to the 
time-averaged Stark shifts disregarded from Hamiltonian ©. 
In principle, the single-ion phases they induce can be compen- 
sated at the end of the evolution by single-ion gates. As for the 
robustness of an implementation however, they are very im- 
portant. Since | Xmi \ is about twenty times as big as 1 mi r] cm | , 
the acquired single-ion phases are much larger than the de- 
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sired conditional n/2 phase, which makes the gate's perfor- 
mance extremely sensitive to laser-intensity fluctuations. One 
way to overcome this is to add additional frequencies onto 
the gate-laser so as to compensate the time-averaged compo- 
nents of the Stark shift without compensating its time- varying 
components i2lll . Alternatively, a spin-echo-technique can be 
applied. In this approach, the gate is divided into two phase- 
space round trips, in each of which half the desired conditional 
phase is acquired, with a pulse interchanging the qubit popu- 
lations in the middle. This increases the gate time by a factor 
of \/2, but for laser power fluctuations slower than the gate 
time both ions pick up the same final single-ion phase. These 
two techniques can also be combined. 

Discussion. Since the scheme is based on the use of weak 
transitions, the large amount of laser power it requires to 
achieve short gate times is somewhat challenging. In a proof- 
of-principle experiment this technical drawback can be cir- 
cumvented using alternative qubit encodings that decrease the 
required laser power. For instance, one can map the qubit onto 
clock states of the hyperfine splitting of the D 5 / 2 manifold. 
Hereby the qubit transition frequency is reduced by almost a 
factor of 500, so that the required laser couplings are reduced 
by more than one order of magnitude. A ground state sub- 
level serves as the Raman mediator and coupling to the other 
Zeeman sublevels of the ground state manifold can be avoided 
by properly choosing the laser polarizations. Since the D$/ 2 
level is now occupied by both qubit states throughout the gate, 
the probability of spontaneous emission is directly given by 
pj = 2 x x T. For instance, for an estimated gate time of 
T = 100 /is (g w 2tt x 285 kHz) including mapping of the 
qubits, the spontaneous scattering per gate-run would still be 
as low as pt = 2 x 10 -4 . Alternatively, another possible en- 



coding is provided by one clock state from S\/2 and another 
from 1^5/2, where values of g w 2n x 210 kHz andpr ~ 10 -4 
are expected i22ll . 

Summary. We propose a new implementation of the (ju- 
gate on hyperfine ground- state qubits that overcomes the two 
main problems of previous implementations: spontaneous 
emission and inefficiency with clock states. The gate is 
driven by forces exerted by two beams in a Raman config- 
uration tuned close to a narrow transition, and state depen- 
dency is achieved by a small frequency detuning instead of 
polarizations. We discuss the idea for 43 Ca + but it can also 
be applied to other ion species with a similar level struc- 
ture. We find a total probability of spontaneous emission per 
gate-run of less than 10 -8 , thus eliminating the main source 
of gate infidelity of previous implementation with magnetic- 
field- sensitive states. In order to achieve short gate times, the 
scheme requires a large amount of laser power. We show two 
examples of how to circumvent this by using alternative qubit 
encodings. With increasingly-powerful laser sources though, 
our (5^ -gate can set a standard for robust high-fidelity entan- 
glement creation in the future large-scale ion-trap quantum 
computer, where the qubit encoding will be that of hyperfine 
ground-state sublevels. With this scheme it seems possible to 
reach the fault tolerant threshold. 
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